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Fluctuations arise universally in Nature as a reflection of the dis- 
crete microscopic world at the macroscopic level. Despite their ap- 
parent noisy origin, fluctuations encode fundamental aspects of the 
physics of the system at hand, crucial to understand irreversibility 
and nonequilibrium behavior. In order to sustain a given fluctuation, 
a system traverses a precise optimal path in phase space. Here we 
show that by demanding invariance of optimal paths under symmetry 
transformations, new and general fluctuation relations valid arbitrar- 
ily far from equilibrium are unveiled. This opens an unexplored route 
toward a deeper understanding of nonequilibrium physics by bring- 
ing symmetry principles to the realm of fluctuations. We illustrate 
this concept studying symmetries of the current distribution out of 
equilibrium. In particular we derive an isometric fluctuation relation 
which links in a strikingly simple manner the probabilities of any 
pair of isometric current fluctuations. This relation, which results 
from the time-reversibility of the dynamics, includes as a particular 
instance the Gallavotti-Cohen fluctuation theorem in this context 
but adds a completely new perspective on the high level of sym- 
metry imposed by time-reversibility on the statistics of nonequilib- 
rium fluctuations. The new symmetry implies remarkable hierarchies 
of equations for the current cumulants and the nonlinear response 
coefficients, going far beyond Onsager's reciprocity relations and 
Green-Kubo formulae. We confirm the validity of the new symme- 
try relation in extensive numerical simulations, and suggest that the 
idea of symmetry in fluctuations as invariance of optimal paths has 
far-reaching consequences in diverse fields. 
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Large fluctuations, though rare, play an important role in 
many fields of science as they crucially determine the fate 
of a system Q]. Examples range from chemical reaction ki- 
netics or the escape of metastable electrons in nanoelectronic 
devices to conformational changes in proteins, mutations in 
DNA, and nucleation events in the primordial universe. Re- 
markably, the statistics of these large fluctuations contains 
deep information on the physics of the system of interest [2] [3] • 
This is particularly important for systems far from equilib- 
rium, where no general theory exists up to date capable of 
predicting macroscopic and fluctuating behavior in terms of 
microscopic physics, in a way similar to equilibrium statisti- 
cal physics. The consensus is that the study of fluctuations 
out of equilibrium may open the door to such general the- 
ory. As most nonequilibrium systems are characterized by 
currents of locally conserved observables, understanding cur- 
rent statistics in terms of microscopic dynamics has become 
one of the main objectives of nonequilibrium statistical physics 
[21 H S E El HIM E M EES El M ED OS QUITS]. Pursuing 
this line of research is both of fundamental as well as practical 
importance. At the theoretical level, the function controlling 
current fluctuations can be identified as the nonequilibrium 
analog of the free energy functional in equilibrium systems 
[21 |3l |H [5] , from which macroscopic properties of a nonequi- 
librium system can be obtained (including its most prominent 
features, as for instance the ubiquitous long range correlations 
[191 120| . etc.) On the other hand, the physics of most mod- 
ern mesoscopic devices is characterized by large fluctuations 
which determine their behavior and function. In this way 
understanding current statistics in these systems is of great 
practical significance. 



Despite the considerable interest and efforts on these is- 
sues, exact and general results valid arbitrarily far from equi- 
librium are still very scarce. The reason is that, while in equi- 
librium phenomena dynamics is irrelevant and the Gibbs dis- 
tribution provides all the necessary information, in nonequi- 
librium physics dynamics plays a dominant role, even in the 
simplest situation of a nonequilibrium steady state [2l [3] |H [5] . 
However, there is a remarkable exception to this absence of 
general results which has triggered an important surge in ac- 
tivity since its formulation in the mid nineties. This is the 
fluctuation theorem, first discussed in the context of simu- 
lations of sheared fluids [15] . and formulated rigorously by 
Gallavotti and Cohen under very general assumptions [16] . 
This theorem, which implies a relation between the probabil- 
ities of a given current fluctuation and the inverse event, is 
a deep statement on the subtle consequences of time-reversal 
symmetry of microscopic dynamics at the macroscopic, irre- 
versible level. Particularly important here is the observation 
that symmetries are reflected at the fluctuating macroscopic 




Fig. 1. The isometric fluctuation relation at a glance. Sketch of the current dis- 
tribution in two dimensions, peaked around its average (J)e, and isometric contour 
lines for different |J|'s. The isometric fluctuation relation, eq. establishes a 

simple relation for the probability of current fluctuations along each of these contour 
lines. 
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level arbitrarily far from equilibrium. Inspired by this illumi- 
nating result, we explore in this paper the behavior of the cur- 
rent distribution under symmetry transformations [21]. Key 
to our analysis is the observation that, in order to facilitate a 
given current fluctuation, the system traverses a well-defined 
optimal path in phase space [3 [3J SI El El El EH • This path is, 
under very general conditions, invariant under certain sym- 
metry transformations on the current. Using this invariance 
we show that for d-dimensional, time-reversible systems de- 
scribed by a locally-conserved field and possibly subject to a 
boundary-induced gradient and an external field E, the prob- 
ability P T (J) of observing a current J averaged over a long 
time r obeys an isometric fluctuation relation (IFR) 
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for any pair of isometric current vectors, |J| = |J'|. Here 
e = e + E is a constant vector directly related to the rate 
of entropy production in the system, which depends on the 
boundary baths via e (see below). 

The above equation, which includes as a particular case 
the Gallavotti-Cohen (GC) result for J' = —J, relates in a 
strikingly simple manner the probability of a given fluctua- 
tion J with the likelihood of any other current fluctuation on 
the d-dimensional hypersphere of radius |J|, see Fig. 1, pro- 
jecting a complex d-dimensional problem onto a much sim- 
pler one-dimensional theory. Unlike the GC relation which is 
a non-differentiable symmetry involving the inversion of the 
current sign, J — s> —J, eq. (TTj) is valid for arbitrary changes 
in orientation of the current vector. This makes the experi- 
mental test of the above relation a feasible problem, as data 
for current fluctuations involving different orientations around 
the average can be gathered with enough statistics to ensure 
experimental accuracy. It is also important to notice that 
the isometric fluctuation relation is valid for arbitrarily large 
fluctuations, i.e. even for the non-Gaussian far tails of cur- 
rent distribution. We confirm here the validity of the new 
symmetry in extensive numerical simulations of two different 
nonequilibrium systems: (i) A simple and very general lattice 
model of energy diffusion 7 8j [23], and (ii) a hard-disk fluid 
in a temperature gradient [24] . 



The Isometric Fluctuation Relation. Our starting point is a 
continuity equation which describes the macroscopic evolution 
of a wide class of systems characterized by a locally-conserved 
magnitude (e.g. energy, particle density, momentum, etc.) 



d t p(r,t) 



Q E [p(r,t)] + £(r,t) 



[2] 



Here p(r,i) is the density field, j(r,t) = QE[p(r,i)] +£(r,t) is 
the fluctuating current, with local average QE[p(r,i)], and 
£(r,t) is a Gaussian white noise characterized by a vari- 
ance (or mobility) a[p(r,t)]. This (conserved) noise term 
accounts for microscopic random fluctuations at the macro- 
scopic level. Notice that the current functional includes in 
general the effect of a conservative external field, Qe[p(i", t)] = 
p(r, f)] + <r[p(r, t)]~E. Examples of systems described by eq. 
2| range from diffusive systems [2E1E1E1E1E1E1E], where 
p(r,t)] is given by Fourier's (or equivalently Fick's) law, 
p(r,t)] = -D[p]Vp(r, f), to most interacting-particle flu- 
ids [251 126] . characterized by a Ginzburg-Landau-type theory 
for the locally-conserved particle density. To completely de- 
fine the problem, the above evolution equation must be sup- 



plemented with appropriate boundary conditions, which may 
include an external gradient. 

We are interested in the probability P T (J) of observing a 
space- and time-averaged empirical current J, defined as 



dt / dr j(r,t). 



T 

This probability obeys a large deviation principle for long 
times [23E8], P T (J) ~ exp[+rL d G(J)], where L is the sys- 
tem linear size and G(J) < is the current large-deviation 
function (LDF), meaning that current fluctuations away from 
the average are exponentially unlikely in time. According to 
hydrodynamic fluctuation theory [21 [4] E] E] , 
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which expresses the locally-Gaussian nature of fluctuations 
EJEIE]- The optimal profile po(r; J) solution of the above vari- 
ational problem can be interpreted as the density profile the 
system adopts to facilitate a current fluctuation J [71 151 [22]. 
To derive eq. (jlj we assumed that (i) the optimal profiles as- 
sociated to a given current fluctuation are time-independent 
[21 El H El El El M El HQ], and (ii) the optimal current 
field has no spatial structure, see Supporting Information (SI). 
This last hypothesis, which greatly simplifies the calculation 
of current statistics, can be however relaxed for our purposes 
(as shown below). The probability Pt(J) is thus simply the 
Gaussian weight associated to the optimal profile. Note how- 
ever that the minimization procedure gives rise to a nonlinear 
problem which results in general in a current distribution with 
non-Gaussian tails [21 El El El El El E] • 

The optimal profile is solution of the following equation 
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stands for functional derivative, and 
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Remarkably, the optimal profile po(r; J) solution of eq. (pH) 
depends exclusively on J and J 2 . Such a simple quadratic de- 
pendence, inherited from the locally-Gaussian nature of fluc- 
tuations, has important consequences at the level of symme- 
tries of the current distribution. In fact, it is clear from eq. 
|5| that the condition 
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implies that po(r;J) will depend exclusively on the magni- 
tude of the current vector, via J 2 , not on its orientation. 
In this way, all isometric current fluctuations characterized 
by a constant |J| will have the same associated optimal pro- 
file, po(r;J) = po(r; |J|), independently of whether the cur- 
rent vector J points along the gradient direction, against it, 
or along any arbitrary direction. In other words, the op- 
timal profile is invariant under current rotations if eq. Q 
holds. It turns out that condition |7| follows from the time- 
reversibility of the dynamics, in the sense that the evolu- 
tion operator in the Fokker-Planck formulation of eq. (|2j 
obeys a local detailed balance condition [1 71 I18j . In this 
case Q E [p(r)]/a[p(r)] = - V5U[p]l5p, with U[p(r)] the sys- 
tem Hamiltonian, and condition (Fn) holds. The invariance of 
the optimal profile can be now used in eq. Q to relate in a 
simple way the current LDF of any pair of isometric current 
fluctuations J and J', with !J| = IJ'I, 



G(J)-G(J') = |e||J|(cos0-cos0') ; 



[81 
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where 9 and 8' are the angles formed by vectors J and J', 
respectively, with a constant vector e = s + E, see below. Eq. 
Q is just an alternative formulation of the isometric fluctu- 
ation relation (fl}. By letting J and J' differ by an infinitesi- 
mal angle, the IFR can be cast in a simple differential form, 
deG(J) = |e||J| sinf?, which reflects the high level of symmetry 
imposed by time-reversibility on the current distribution. 

The condition 8wi[p(r)]/Sp(r') = can be seen as a con- 
servation law. It implies that the observable wi[p(r)] is in 
fact a constant of motion, e = wi[/s(r)], independent of the 
profile p(r), which can be related with the rate of entropy 
production via the Gallavotti-Cohen theorem |16l 1171 ITS] . In 
a way similar to Noether's theorem, the conservation law for 
e implies a symmetry for the optimal profiles under rotations 
of the current and a fluctuation relation for the current LDF. 
This constant can be easily computed under very general as- 
sumptions (see SI). 

Implications and Generalizations. The isometric fluctuation 
relation, eq. |TJ, has far-reaching and nontrivial consequences. 
First, the IFR implies a remarkable hierarchy of equations for 
the cumulants of the current distribution, see eq. ( |13[ ) in 
Methods. This hierarchy can be derived starting from the 
Legendre transform of the current LDF, p(A) = maxj[G(J) + 
A • J] , from which all cumulants can be obtained [3] , and writ- 
ing the IFR for //(A) in the limit of infinitesimal rotations. As 
an example, the cumulant hierarchy in two dimensions implies 
the following relations 

{J x } e = TL 2 [e x (AJ 2 }^e y (AJ x AJ y ) t ] [9] 
(J y ) e = TL 2 [e y {AJ 2 ).-e x {AJ x AJ y } c ] 



[101 



2(AJ x AJ y ) c = TL 2 [e y {Ajl) t -e x (AJ 2 x AJy) e 
= tL 2 [ Ei (AJ s 3 ) £ - e y (AJ x AJ 2 ) c ] 
(AJ^)e - <AJ y 2 ) E = tL 2 [e x (AJ x AJ 2 } ( - e y (AJ*AJ v ) e ] , 

for the first cumulants, with AJ a = J a — (J a )e- It is 
worth stressing that the cumulant hierarchy is valid arbitrar- 
ily far from equilibrium. In a similar way, the IFR implies a 
set of hie r arch ies for the nonlinear response coefficients, see 
eqs. ( |15[ )-(|17| in Methods. In our two-dimensional exam- 

Pie, let {n) X ( „ x , nH ) 

lant {AJ X * AJy "} e to order e x x e v v , witn n — n x t il v 
k — k x + k y . To the lowest order these hierarchies imply On- 
sager's reciprocity symmetries and Green-Kubo relations for 
the linear response coefficients of the current. They further 
predict that in fact the linear response matrix is proportional 
to the identity, so § x ™ = $ x [^> - <°>~<°.°> - <°^°> 

while w x gg = = w x g;;j. 

of the current can be simply written in terms of the linear 
coefficients of the second cumulants as 



be the response coefficient of the cumu- 
AJy v } c to order e x x ey v , with n — n x + n y and 
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(i)X(i;J) = 2 [(2)X(°;J) + (2)X(?;i)] (symmetric results hold for 
n x = 0, n y = 1). Linear response coefficients for the second- 
order cumulants also obey simple relations, e.g. pjXn'i) = 

(i)„(o,i) , (1) (1,0) (1) (0,1) _ (1) (1,0) (1) (0,1)' , 
dllQ ( 2 )X(2,0) + (2)X(2,0) — (2)X(o,2) + (2)A(o,2)' dI1U 



(2)X(i,i) 

the set of relations continues to arbitrary high orders. In this 
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way hierarchies ( 15 l-( 17 1 , which derive from microreversibility 
the IFR, provide deep insights into nonlinear 



response theory for nonequilibrium systems [29] . 

The IFR and the above hierarchies all follow from the 
invariance of optimal profiles under certain transformations. 
This idea can be further exploited in more general settings. 



In fact, by writing explicitly the dependence on the external 
field E in eq. <f5t for the optimal profile, one realizes that 
if j^rj J Q[p( r )]dr = 0, together with the time-reversibility 
condition, eq. ffb, the resulting optimal profiles are invariant 
under independent rotations of the current and the external 
field. It thus follows that the current LDFs for pairs (J, E) 
and (J' = 1ZJ,E* = SE), with 1Z, S independent rotations, 
obey a generalized isometric fluctuation relation 

C?e(J)-Ge*(J') = e-(J-j')-^-(E-E*)+J-E-j'-E* , [11] 

where we write explicitly the dependence of the current LDF 
on the external field. The vector v = J Q[p(r)]dr is now an- 
other constant of motion, independent of p(r), which can be 
easily computed (see SI). For a fixed boundary gradient, the 
above equation relates any current fluctuation J in the pres- 
ence of an external field E with any other isometric current 
fluctuation J' in the presence of an arbitrarily-rotated exter- 
nal field E*, and reduces to the standard IFR for E = E*. 
Condition -g-jzr\ J Q{p{ r )]dr = is rather general, as most 
time-reversible systems with a local mobility o~[p] do fulfill 
this condition (e.g., diffusive systems). 

The IFR can be further generalized to cases where the 
current profile is not constant, relaxing hypothesis (ii) above. 
Let P T [^7"(r)] be the probability of observing a time-averaged 
current field «7(r) = r _1 J Q T dt j(r, t). This vector field must 
have zero divergence because it is coupled via the continuity 
equation to an optimal density profile which is assumed to be 
time-independent, see SI and hypothesis (?) above. Because 
of time-reversibility, QE[p(r)]/cr[p(r)] = — W5H[p]/Sp and it 
is easy to show in the equation for the optimal density profile 
that the term linear in ^T(r) vanishes, so po[v; remains 
invariant under (local or global) rotations of J7"(r), see SI. In 
this way, for any divergence- free current field «7'(r) locally- 
isometric to «7"(r), so ^7"'(r) 2 = 3(?) 2 Vr, we can write a 
generalized isometric fluctuation relation 



lim 



hi 



Pr[J-(r) 



Pr[^'(r)] 



[12] 

where the integral (whose result is independent of p(r)) is 
taken over the boundary dA of the domain A where the sys- 
tem is defined, and h is the unit vector normal to the bound- 
ary at each point. Eq. ( 12 1 generalizes the IFR to situations 



where hypothesis (ii) is violated, opening the door to isome- 
tries based on local (in addition to global) rotations. As a 
corollary, we show in the SI appendix that a similar general- 
ization of the isometric fluctuation symmetry does not exist 
whenever optimal profiles become time-dependent, so the IFR 
breaks down in the regime where hypothesis (i) is violated. 
In this way, we may use violations of the IFR and its gen- 
eralizations to detect the instabilities which characterize the 
fluctuating behavior of the system at hand [2l l9l 110] . 

Checking the Isometric Fluctuation Relation. We have tested 
the validity of the IFR in extensive numerical simulations of 
two different nonequilibrium systems. The first one is a sim- 
ple and very general model of energy diffusion [7| [H] de- 
fined on a two-dimensional (2D) square lattice with L 2 sites. 
Each site is characterized by an energy a, i S [1,L 2 ], and 
models a harmonic oscillator which is mechanically uncou- 
pled from its nearest neighbors but interact with them via a 
stochastic energy-redistribution process. Dynamics thus pro- 
ceeds through random energy exchanges between randomly- 
chosen nearest neighbors. In addition, left and right boundary 
sites may interchange energy with boundary baths at tempera- 
tures Tl and Tr, respectively, while periodic boundary condi- 
tions hold in the vertical direction. For Tl 7^ Tr the systems 
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reaches a nonequilibrium steady state characterized, in the 
absence of external field (the case studied here), by a linear 
energy profile p s t(r) — Tl+x (Tr — Tl) and a nonzero average 
current given by Fourier's law. This model plays a fundamen- 
tal role in nonequilibrium statistical physics as a testbed to 
assess new theoretical advances, and represents at a coarse- 
grained level a large class of diffusive systems of technological 
and theoretical interest [7] [HI. The model is described at the 
macroscopic level by eq. pj with a diffusive current term 
Q[p(r,t)] = -D[p]Vp with~D[p] = | and a[p] = p 2 , and it 
turns out to be an optimal candidate to test the IFR because: 
(1) the associated hydrodynamic fluctuation theory can be 
solved analytically [30], and (2) its dynamics is simple enough 
to allow for a detailed numerical study of current fluctuations. 

In order to test the IFR in this model we performed a large 
number of steady-state simulations of long duration r > L 2 
(the unit of time is the Monte Carlo step) for L = 20, Tl — 2 
and Tr = 1, accumulating statistics for the space- and time- 
averaged current vector J. The measured current distribution 
is shown in the bottom inset to Fig. 2, together with a fine 
polar binning which allows us to compare the probabilities 
of isometric current fluctuations along each polar corona, see 
eq. Q. Taking G(J) = (rL d ) _1 lnP r (J), Fig. 2 confirms 
the IFR prediction that G(J) — G(J'), once scaled by J| _1 , 
collapses onto a linear function of cos — cos 0' for all values 
of |J|, see eq. (JsJ) - Here 8, 6' are the angles formed by the 
isometric current vectors J, J' with the i-axis (E = in our 
case) . We also measured the average energy profile associated 
to each current fluctuation, po(r;«J), see top inset to Fig. 2. 
As predicted above, profiles for different but isometric cur- 
rent fluctuations all collapse onto a single curve, confirming 
the invariance of optimal profiles under current rotations. 

Standard simulations allow us to explore moderate fluctu- 
ations of the current around the average. In order to test the 
IFR in the far tails of the current distribution, corresponding 
to exponentially unlikely rare events, we implemented an el- 
egant method recently introduced to measure large deviation 
functions in many-particle systems [32]. The method, which 
yields the Legendre transform of the current LDF, p(X), is 
based on a modification of the dynamics so that the rare 
events responsible of the large deviation are no longer rare 
[32], and has been recently used with success to confirm an 
additivity conjecture regarding large fluctuations in nonequi- 
librium systems [Tj [8] . Using this method we measured p.(X) 
in increasing manifolds of constant |A+e|, see Fig. 3. The IFR 
implies that p(X) is constant along each of these manifolds, 
or equivalently p(A) = fj,[R,$(X + e) — e], V0 £ [0, 2tt], with 
72,0 a rotation in 2D of angle <f>. Fig. 3 shows the measured 
/x(A) for different values of |A + e| corresponding to very large 
current fluctuations, different rotation angles (j> and increasing 
system sizes, together with the theoretical predictions [3D] . As 
a result of the finite, discrete character of the lattice system 
studied here, we observe weak violations of IFR in the far tails 
of the current distribution, specially for currents orthogonal to 
e. These weak violations are expected since a prerequisite for 
the IFR to hold is the existence of a macroscopic limit, i.e. eq. 
|2| should hold strictly, which is not the case for the relatively 
small values of L studied here. However, as L increases, a clear 
convergence toward the IFR prediction is observed as the ef- 
fects associated to the underlying lattice fade away, strongly 
supporting the validity of IFR in the macroscopic limit. 

We also measured current fluctuations in a Hamiltonian 
hard-disk fluid subject to a temperature gradient ,24 . This 
model is a paradigm in liquid state theory, condensed matter 
and statistical physics, and has been widely studied during 
last decades. The model consists in N hard disks of unit di- 
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Fig. 2. Confirmation of IFR in a diffusive system. The IFR predicts that 
| J | — 1 [G( J) — G(J')] collapses onto a linear function of COS 6 — COS 0' for 
all values of |J|. This collapse is confirmed here in the energy diffusion model for a 
wide range of values for |J|. Bottom inset: Measured current distribution together 
with the polar binning used to test the IFR. Top inset: Average profiles for different 
but isometric current fluctuations all collapse onto single curves, confirming the in- 
variance of optimal profiles under current rotations. Angle range is \6\ < 16.6°, 
see marked region in the histogram. 
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Fig. 3. IFR for large current fluctuations. Legendre transform of the cur- 
rent LDF for the energy diffusion model, for different values of |A + e| corre- 
sponding to very large current fluctuations, different rotation angles such that 
A' = (A + e) — e, and increasing system sizes. Lines are theoretical pre- 
dictions. The IFR predicts that fl(X) = /i[7^(A + e) - e] V0 6 [0, 2tt]. 
The isometric fluctuation symmetry emerges in the macroscopic limit as the effects 
associated to the underlying lattice fade away. 



ameter interacting via instantaneous collisions and confined 
to a box of linear size L such that the particle density is fixed 
to $ = N/L 2 = 0.58. Here we choose N = 320. The box is 
divided in three parts: a central, bulk region of width L — 2a 
with periodic boundary conditions in the vertical direction, 
and two lateral stripes of width a — L/4 which act as de- 
terministic heat baths, see bottom inset to Fig. 4. This is 
achieved by keeping constant the total kinetic energy within 
each lateral band via a global, instantaneous rescaling of the 



1 Notice that p-dependent corrections to a constant mobility a , which are typically irrelevant from 
a renormalization-group point of view 26 , turn out to be essential for current fluctuations as they 
give rise to non-Gaussian tails in the current distribution. 
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Fig. 4. IFR in a hydrodynamic hard-disk fluid. Confirmation of IFR in a two di- 
mensional hard-disk fluid under a temperature gradient after a polar binning of the 
measured current distribution. As predicted by IFR, the difference of current LDFs 
for different isometric current fluctuations, once scaled by the current norm, collapses 
in a line when plotted against COS — cosf?'. Top inset: Optimal temperature 
profiles associated to different current fluctuations. Profiles for a given |J| and dif- 
ferent angles 9 £ [ — 7.5°, +7.5°] all collapse onto a single curve, thus confirming 
the invariance of optimal profiles under current rotations. Notice that the profiles 
smoothly penetrate into the heat baths. Bottom inset: Snapshot of the 2D hard-disk 
fluid with Gaussian heat baths. 



velocity of bath particles after bath-bulk particle collisions. 
This heat bath mechanism has been shown to efficiently ther- 
mostat the fluid [24], and has the important advantage of 
being deterministic. As for the previous diffusive model, we 
performed a large number of steady state simulations of long 
duration (r > 2N collisions per particle) for Tl = 4 and 
Tr — 1, accumulating statistics for the current J and mea- 
suring the average temperature profile associated to each J. 
Fig. 4 shows the linear collapse of |J| _1 [G(J) — G(J')] as a 
function of cosf? — cosf?' for different values of |J|, confirming 
the validity of the IFR for this hard-disk fluid in the moderate 
range of current fluctuations that we could access. Moreover, 
the measured optimal profiles for different isometric current 
fluctuations all nicely collapse onto single curves, see top inset 
to Fig. 4, confirming their rotational invariance. 

It is interesting to notice that the hard-disk fluid is a fully 
hydrodynamic system, with 4 different locally-conserved cou- 
pled fields possibly subject to memory effects, defining a far 
more complex situation than the one studied here, see eq. q2p . 
Therefore the validity of IFR in this context suggests that tnis 
fluctuation relation, based on the invariance of optimal profiles 
under symmetry transformations, is in fact a rather general 
result valid for arbitrary fluctuating hydrodynamic systems. 

A few remarks are now in order. First, as a corollary 
to the IFR, it should be noted that for time-reversible sys- 
tems with additive fluctuations, i.e. with a constant, profile- 
independent mobility a, the optimal profile associated to a 

§en current fluctuation is in fact independent of J, see eq. 
, and hence equal to the stationary profile. In this case it 
is easy to show that current fluctuations are Gaussian, with 
G(J) = e • (J — (J) e ) + rr-^J 2 - (J) 2 ). This is the case, for 
instance, of model B in the Hohenberg-Halperin classification 
|26| Q On the other hand, it should be noticed that the time- 
reversibility condition for the IFR to hold, eq. |7| , is just 
a sufficient but not necessary condition. In fact, we cannot 
discard the possibility of time-irreversible systems such that, 
only for the optimal profiles, SLJi[p(r)]/Sp(r')\ P0 = 0. 



Discussion. The IFR is a consequence of time-reversibility for 
systems in the hydrodynamic scaling limit, and reveals an un- 
expected high level of symmetry in the statistics of nonequilib- 
rium fluctuations. It generalizes and comprises the Gallavotti- 
Cohen fluctuation theorem for currents, relating the proba- 
bilities of an event not only with its time-reversal but with 
any other isometric fluctuation. This has important conse- 
quences in the form of hierarchies for the current cumulants 
and the linear and nonlinear response coefficients, which hold 
arbitrarily far from equilibrium and can be readily tested in 
experiments. A natural question thus concerns the level of 
generality of the isometric fluctuation relation. In this paper 
we have demonstrated the IFR for a broad class of systems 
characterized at the macroscale by a single conserved field, us- 
ing the tools of hydrodynamic fluctuation theory (HFT). This 
theoretical framework, summarized in the path large devia- 
tion functional, eq. (3) in the SI appendix, has been rigorously 
proven for a number of interacting particle systems [21 13 SI E] t 
but it is believed to remain valid for a much larger class of 
systems. The key is that the Gaussian nature of local fluc- 
tuations, which lies at the heart of the approach, is expected 
to emerge for most situations in the appropriate macroscopic 
limit as a result of a central limit theorem: although micro- 
scopic interactions can be extremely complicated, the ensu- 
ing fluctuations of the slow hydrodynamic fields result from 
the sum of an enormous amount of random events at the mi- 
croscale which give rise to Gaussian statistics. There exist 
of course anomalous systems for which local fluctuations at 
the macroscale can be non-Gaussian. In these cases we can- 
not discard that a modified version of the IFR could remain 
valid, though the analysis would be certainly more compli- 
cated. Furthermore, our numerical results show that the IFR 
remains true even in cases where it is not clear whether the 
HFT applies, strongly supporting the validity of this symme- 
try for arbitrary fluctuating hydrodynamic systems. 

A related question is the demonstration of the IFR start- 
ing from microscopic dynamics. Techniques similar to those 
in Refs. |17l I31j . which derive the Gallavotti-Cohen fluctua- 
tion theorem from the spectral properties of the microscopic 
stochastic evolution operator, can prove useful for this task. 
However, in order to prove the IFR these techniques must 
be supplemented with additional insights on the asymptotic 
properties of the microscopic transition rates as the macro- 
scopic limit is approached. In this way we expect finite-size 
corrections to the IFR which decay with the system size, as it 
is in fact observed in our simulations for the energy diffusion 
model, see Fig. 3. Also interesting is the possibility of an 
IFR for discrete isometries related with the underlying lattice 
in stochastic models. These open questions call for further 
study. 

We have shown in this paper how symmetry principles 
come forth in fluctuations far from equilibrium. By demand- 
ing invariance of the optimal path responsible of a given 
fluctuation under symmetry transformations, we unveiled a 
novel and very general isometric fluctuation relation for time- 
reversible systems which relates in a simple manner the proba- 
bility of any pair of isometric current fluctuations. Invariance 
principles of this kind can be applied with great generality 
in diverse fields where fluctuations play a fundamental role, 
opening the door to further exact and general results valid ar- 
bitrarily far from equilibrium. This is particularly relevant in 
mesoscopic biophysical systems, where relations similar to the 
isometric fluctuation relation might be used to efficiently mea- 
sure free-energy differences in terms of work distributions [33] . 
Other interesting issues concern the study of general fluctua- 
tion relations emerging from the invariance of optimal paths 
in full hydrodynamical systems with several conserved fields, 
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or the quantum analog of the isometric fluctuation relation in 
full counting statistics. 

Appendix: Hierarchies for the cumulants and response co- 
efficients 

The moment-generating function associated to P r (J), de- 
fined as II T (A) = J P T (J) exp(rL d A • J)dJ, scales for long 
times as n T (A) ~ exp[+rL d /u(A)], where /i(A) = maxj[G(J) + 
A ■ J] is the Legendre transform of the current LDF. The 
cumulants of the current distribution can be obtained from 
the derivatives of /Ct(A) evaluated at A = 0, i.e. p|™' , = 

[9>(A)W...A^] M = (TL <i ) n - 1 <A^ 1 ...AJj*) e 1 far n > 
1, where AJ a = J a — (1 — S ni i){J a )e and <5 n , m is the Kronecker 
symbol. The IFR can be stated for the Legendre transform of 
the current LDF as fl(X) — p[7Z(\ + e) — e], where 1Z is any 
d-dimensional rotation. Using this relation in the definition of 
the n-th order cumulant in the limit of infinitesimal rotations, 
TZ = I + A9£, it is easy to show that 



naL -P°>L L (n 1 ...n a -l...n f3 + l...n d ) + e " i 'T'M(n 



(n+1) 



■«d) 



o, 

[13] 

where £ is any generator of d-dimensional rotations, and sum- 
mation over repeated Greek indices (G [1, d]) is assumed. The 
above hierarchy relates in a simple way cumulants of orders n 
and n+1 Vn > 1, and is vali d arb itrarily far from equilibrium. 
As an example, eqs. fl9p and ([Top above show the first two sets 
of relations (n = 1,2) of the above hierarchy in two dimen- 
sions. In a similar way, we can explore the consequences of 
the IFR on the linear and nonlinear response coefficients. For 
that, we now expand the cumulants of the current in powers 
of e 



{n 1 ...n d ) V 



oo ^ k 



(.k) lkx...k d ) fci fe d 
fcl (n)X(m... ni ) e i - e d 

fc=0 fe 1 ...fe d =0 

Ei k t =k 



[14] 



Inserting expansion ( 14 1 into the cumulant hierarchy, eq. ( 13 1, 
and matching order by order in k, we derive another inter- 
esting hierarchy for the response coefficients of the different 
cumulants. For k = this reads 



r (0) (0...0) 

n a L,p a ( 7l )X(„ 1 ... no 



-l...na+l...n d ) 



0. 



[151 



which is a symmetry relation for the equilibrium (e = 0) cur- 
rent cumulants. For k > 1 we obtain 



E 



ki...k d =0 

J2i fci=fc>l 



' n (k) (kx ■■■k d ) 
t ~P a (n)X-(nx...n a - 



X...np+l...n d ) 



j.r (k-l) (k 1 ...k„-l...k i ) 

Tinrw („+l)A( ni ... nT+ l... n<J ) 



= 0, 



[161 



which relates fc-order response coefficients of n-order cumu- 
lants with (k — l)-order coefficients of (n + l)-order cumu- 
lants. Relations ( 15 )-( 16 ) for the response coefficients re- 
sult from the IFR m the limit of infinitesimal rotations. For 
a finite rotation 1Z = —I, which is equivalent to a current 
inversion, we have p(A) = p(— A — 2e) and we may use 

this in the definition of response coefficients, ^ Xfa' nti = 

k\{d n+k f i(\)/d\ , l 1 ...\yde k 1 1 ...e k d d } x = 0= e, see eq. Q, to ob- 
tain a complementary relation for the response coefficients 



fei k d 

(fe) (kx...k d ) , , v-^ 

(n)X( ni ...n d ) Z^t (k-v) 

P1 =0 p d =0 yJ 



{-l) n+P 2 P (fc_ p ) (k 1 -p 1 ...k d -p d ) 

(n+p)^-(n 1 +p 1 ...n d +p d ) ' 

[17] 



where p = Y2iPi' A similar equation was derived in [29] 
from the standard fluctuation theorem, a lthough the IFR adds 
further relations. All together, eqs. (15l-(17l imply deep 



relations between the response coefficients at arbitrary or- 
ders which go far beyond Onsager's reciprocity relations and 
Green-Kubo formulae. As an example, we discuss in the main 
text some of these relations for a two-dimensional system. 

Appendix: Hydrodynamic fluctuation theory 

The evolution of the system of interest is described by the 
following Langevin equation 



d t p(r,t) 



Q E [p(r,t)]+€(r,t) 



[181 



which expresses the local conservation of certain physi- 
cal observable. Here p(r, t) is the density field, j(r, t) = 
QE[p(r,t)] + £(r,i) is the fluctuating current, with local av- 
erage QE[p(r, i)], and £(r,i) is a Gaussian white noise with 
zero mean and characterized by a variance (or mobility) 
a[p(r,t)]. Notice that the current functional includes in gen- 
eral the effect of a conservative external field, QE[p(r,f)] = 
Q[p(r, t)] + rr[p(r, i)]E. Using a path integral formulation [30] . 
the probability of observing a given history {p(r, t), j(r, t)}g 
of duration t for the density and current fields can be written 



P({pJ}S)~exp( + L d / T [p,j] 



[191 



where L is the system linear size, d is the dimensionality, and 
the functional I T [p, j] is 



It [pj] 



dt / dv 



j(r,t)-Q E [p(r,t); 



201 



, j 2a[p(r,t)] 

with p(r, t) and j(r, t) coupled via the continuity equation 

r3 t p(r,t) + V-j(r,t)=0. [21] 

In this way the probability of each history {p,j}o has a 
Gaussian weight around the aver age local behavior given by 
Qe[p(i", t)\. Eqs. (191 and (20 1 are equivalent to the hy- 



drodynamic fluctuation theory recently proposed by Bertini 
and coworkers [U O [2] . The probability P T ( J) of observing a 
space- and time-averaged empirical current J, defined as 



dt I dr j(r,t) . 



[221 



can be obtained from the path integral of P ({p, j}o) restricted 
to histories {p, j}o compatible with a given J, 



VpVjP({ P ,i} T )s(j 



dt / dr j(r,t) 



P T (J) : 

J V T Jo J 

[23] 

This probability scales for long times as P T (J) ~ 
exp[+rL d G(3)], and the current large deviation function 
(LDF) G(J) can be related to I T [p, j] via a simple saddle-point 
calculation in the long-time limit, 



G(J) = hm — max 7 r [p,j] 

T— >CX> T P(r,t) 
j(r.t) 



[241 



subject to constraints (21 1 and (22 1. The density and cur 



rent fields solution of this variational problem, denoted here 
as po(r, t; J) and jo(r, t; J), can be interpreted as the optimal 
path the system follows in order to sustain a lont-time cur- 
rent fluctuation J. It is worth emphasizing here that the ex- 
istence of an optimal path rests on the presence of a selection 
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principle at play, namely a long time, large size limit which 
selects, among all possible paths compatible with a given fluc- 
tuation, an optimal one via a saddle point mechanism. Eq. 
( |24[ | defines a complex spatiotemporal problem whose solution 
remains challenging in most cases [30l l4l f5l f2l l6l fTl l8l l9l flO] . 
However, the following hypotheses greatly reduce its complex- 
ity: 

(i) We assume that the optimal profiles responsible of a given 
current fluctuation are time-independent [6], po(r; J) and 
jo(r;J). This, together with the continuity equation, im- 
plies that the optimal current vector field is divergence-free, 
V-j (p;J) = 0. 

(ii) A further simplification consists in assuming that the op- 
timal current field is in fact constant across space, so 
jo(r; J) = J. 

Provided that these hypotheses hold, the current LDF can be 
written as 



G(J) = - min 



(J - Q E [p(r)]r 



P (r) J 2a[p(r) 



dr . 



[251 



The optimal density profile is thus solution of the following 
differential equation 

<Wp( r )l 2 j <Wl(r)] | j2 <Wp(r)] = Q , 2Q , 



6p(r>) 



6p(r>) 



6p(r>) 



where , , ,, stands for functional derivative, and 

u, n [p(r)] = / drW„[p(r)] with W«[p(r)] = Q ^ ( ^ . 

[27] 

For time-reversible systems, one can see that the ev olution 
operator in the Fokker-Planck formulation of eq. ( 18 1 obeys 
a local detailed balance condition, and 



W 1 [p(r)] = 



Q E [p(r) 
a[p(r)] 



8p 



[281 



where H[p(r)] is the system Hamiltonian. In this case, by 
using vector integration by parts, it is easy to show that 



dF^^(r).n = 0, 
bp 

[291 



for any divergence- free vector field «4.(r). The second inte- 
gral is taken over the boundary <9A of the domain A where 
the system is defined, and h is the unit vector normal to the 
boundary at e ach point. In particular, by taking «4.(r) = J 
constant, eq. ( 29 1 implies that 8wi[p(r)]/5p(r') = 0. Hence 
for time-reversible systems the optimal profile po(r; J ) rem ains 
invariant under rotations of the current J, see eq. (261, and 



this allows us to prove the isometric fluctuation relation (IFR), 
eqs. (1) and (8) in the main text. 

We can now relax hypothesis (ii) above and study cases 
where the current profile is not constant. Let P r [^7"(r)] be 
the probability of observing a time-averaged current field 
«7(r) = t -1 J Q T dt j(r, t). Notice that this vector field must 
be divergence-free because of hypothesis (i). This proba- 
bility also obeys a large deviation principle, P T [j7"(r)] ~ 
cxp (+rL d G[j7"(r)]) , with a current LDF equivalent to that 
in eq. ( |25[ ) but with a space-dependent current field 
The optimal density profile po[r; ^"( r )] is now solution of 



5p(r>) 



dr (lV 2 [p(r)]-2J'(r)-W 1 [ /9 (r)]+J' 2 (r)VKo[p(r)]) =0, 



which is the equivalent to eq. (|26[ ) in this case. For time- 
reversible systems condition |2ij| holds and po[r;j7"(r)] re- 
mains invariant under (local or global) rotations of Sir). In 
this way we can simply relate P T [j7'(r)] with the probability of 
any other divergence-free current field «7"'(r) locally-isometric 
to J(r), i.e. J\rf = J{r) 2 Vr, via a generalized isometric 
fluctuation relation, see eq. (12) in the paper. Notice that in 
general an arbitrary local or global rotation of a divergence- 
free vector field does not conserve the zero-divergence prop- 
erty, so this constraints the current fields and/or local rota- 
tions for which this generalized IFR applies. 

The large deviation function for the space- and time- 
averaged current, G(J), can be related to G[S(r)] via a con- 
traction principle 



G(J) = max G[J(r)] . 

.7(r):V J-(r) = 



[311 



The optimal, divergence-free current field J (r;J) solution of 
this variational problem may have spatial structure in general. 
However, numerical results and phenomenological arguments 
strongly suggest that the constant solution, j7"o(r; J) = J, is 
the optimizer at least for a wide interval of current fluctua- 
tions, showing that hypothesis (ii) above is not only plausible 
but also well justified on physical grounds. In any case, the 
range of validity of this hypothesis can be explored by study- 
ing the limit of local stability of the constant current solution 
using tools similar to those in Ref. [9]. 

Hypotheses (i) and (ii) are the straightforward general- 
ization to d-dimensional systems of the Additivity Princi- 
ple recently conjectured by Bodineau and Derrida for one- 
dimensional diffusive systems [6]. This conjecture, which has 
been recently confirmed for a broad current interval in exten- 
sive simulations of a general diffusion model 7, 8 , is however 
known to break down in some special cases for extreme cur- 
rent fluctuations, where time-dependent profiles in the form of 
traveling waves propagating along the current direction may 
emerge 1 5 2 !) ID . As in previous cases, we can now study 
the probability P ({j(r, t)}5) of observing a particular history 
for the current field, whi ch c an be written as the path integral 
of the probability in eq. ( 19 1 over histories of the density field 
{p(r, t) }g c oupled to the desired current field via the continu- 
ity eq. ( 21 1 at every point in space and time. This probability 



obeys another large deviation principle, with an optimal his- 
tory of the density field {po (r, t)}p which is solution of an 
equation similar to eq. ( 30 1 but with time-dependent pro- 
files. However, as opposed to the cases above, the current field 
j(r, t) is not necessarily divergence-free because of the time- 
dependence of th e associated po(r,t), resulting in a violation 
of condition ( 29 1 . In this way the optimal po (r, t) depends on 
both j(r, t) and j(r, t) 2 so it does not remain invariant under 
(local or global) instantaneous rotations of the current field, 
resulting in a violation of the generalized isometric fluctuation 
relation in the time-dependent regime. 

Notice that the dynamic phase transition to time- 
dependent optimal paths is expected to occur only for extreme 
current fluctuations, thus rendering valid the isometric fluc- 
tuation relations for a wide, subcritical current interval. In- 
terestingly, we can use the IFR to detect such dynamic phase 
tran sitio n. If we measure P T [^7"(r)] in a system described by 
eq. ( 18 1 at the macroscopic level, finding that the measured 



[301 



probabilities do not obey the generalized IFR, then we can 
conclude that such a violation of IFR is due to the onset of 
time-dependent optimal profiles, thus signaling the dynamic 
phase transition. On the other hand, breakdown of the stan- 
dard IFR (for space- and time-averaged currents) may signal 
the onset of space-dependent, divergence-free optimal current 
profiles or the aforementioned dynamic phase transition. In 



Footline Author PNAS | Issue Date | Volume | Issue Number | 7 



this way, the combined use of the IFR and its generalizations 
is capable of a full characterization of the instabilities which 
characterize the fluctuating behavior of the system at hand 

Gaum. 



[321 



Appendix: Constants of motion 

A sufficient condition for the IFR to hold is that 
<5wi[p(r)] 



Sp(r>) 



= 0, 



with the functional ui[p(r)] d efine d in eq. (27 1 above. 
We have shown that condition ( 32 1 follows from the time- 
reversibility of the dynamics, in the sense that the evolution 
operator in the Fokker-Planck formulati on o f eq. pf obe ys a 
local detailed balance condition, see eq. (291. Condition (32 l 
implies that wi[p(r)] is in fact a constant of motion, e, inde- 
pendent of the profile p(r). Therefore we can use an arbitrary 
profile p(r), compatible with boundary conditions, to compute 
e. We now choose boundary conditions to be gradient-like in 
the ^-direction, with densities p.L and pn at the left and right 
reservoirs, respectively, and periodic boundary conditions in 
all other directions. Given these boundaries, we now select a 
linear profile 

p(r) = p L + (pr - Pl)x, [33] 

to compute e, with x € [0,1], and assume very general forms 
for the current and mobility functionals 

Q[p(r)] = A,,o[p]V P + A lm [p](V m p) 2 "Vp, 

n,m>0 

cr[p(r)] = cr , [p]+ Yl a n m [p](^ m pf n , 

n,m>0 

where as a convention we denote as F[p] a generic functional 
of the profile but not of its derivatives. It is now easy to show 
that e — ex + E, with 



g= f PR d DoAp) + E„>o DrajpKpR - pl) 2 

J PL °"0,o(p) + Em>0 °" ml 



[341 



and x the unit vector along the gradient direction. In a similar 
way, if the following condition holds 



8 

Mr 7 ) 



Q[p(r)]dr = 0, 



[351 



together with time-reversibility, eq. (32 1, the system can 
be shown to obey an extended isometric fluctuation relation 
which links any current fluctuation J in the presence of an 
external field E with any other isometric current fluctuation 
J' in the presence of an arbitrarily-rotated external field E* , 
and reduces to the stand ard IFR for E = E*, see eq. (11) 
in the paper. Condition (351 implies that v = J Q[p(r)]dr 
is another constant of motion, which can be now written as 
v = vx, with 



dp 



D ,o(p) + Y, d ^(p)(PR 



PL) 



[361 



As an example, for a diffusive system Q[p(r)] = — D[p]Vp(r), 
with D[p] the diffusivity functional, and the above equations 
yield the familiar results 



i-pl D 

J PR ° 



D(p± 
(p) 



dp, 



D(p)dp, 



for a standard local mobility o[p]. 
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